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The structural integrity of space flight hardware is established by 
a combination of qualification tests and analyses which simulate actual 
operating conditions, including flight loads, temperatures, and corrosive 
environments. These structural analysis and test activities usually 
fall into three distinct areas. The first two areas, strength and 
fatigue analysis, assume the load carrying structure is unflawed. This 
assumption implies that no defects have been introduced during the 
manufacturing process of each individual part which in reality, can 
never be possible on an economical basis. 

The existence of flaws is accounted for in the third area, fracture 
mechanics. This area becomes an important effort in which defects are 
known as a result of quality inspections, or assumed to exist in a part 
and an assessment is made as to their impact on the parts useful life. 
Fracture mechanics attempts to predict the useful service life of an 
initially flawed structural part by calculating crack growth and eventual 
part failure due to unstable crack growth. 

This paper compares the service life calculations of two computer 
codes, NASCRAC and NASA/FLAGRO. The analysis technique is based on 
linear elastic fracture mechanics (LEFM) , in which stresses remain 
below the yield strength of an elastic/plastic material. Subcritical 
crack growth calculations assume that in a metallic part, the extent of 
yielding at the crack tip is very small compared to the crack size and 
the bulk of the cracked body remains elastic. 

To perform service life calculations, one must have a relationship 
expressing incremental crack growth, DA/DN, as a function of loading, 
geometry, and material properties. Load and geometry are expressed in 
terms of the cyclic stress intensity factor, AK. The crack growth rate 
as a function of AK is then determined by material tests, plotting DA/DN 
versus AK for the given material, loading condition, and environment. 

Crack growth rate equations such as the Paris, Walker, and modified 
Forman equations are used to obtain a “best fit” curve to the laboratory 
DA/DN versus AK data. Constants in the equations which result in a 
“best fit" then become crack growth rate material constants for a par- 
ticular set of laboratory conditions. 

Two extreme values of AK also become material constants; Is 

the threshold stress intensity, below which no crack growth occurs, K^ 
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is the critical stress intensity at which a crack becomes unstable and 
complete fracture occurs. Formulations of AK solutions and crack growth 
rate equations form the basis of computer codes which numerically inte- 
grated the DA/DN=*F(AK) relationship. 

Before a computer code is used as part of the structural integrity 
assessment process, it should be exercised thoroughly and its numerical 
calculations checked to insure reasonable and accurate answers. The 
results presented herein compare the safelife calculations of two com- 
puter codes with each other, and with test data to a limited extent. 

The computer program, NASA/FLAGRO (commonly known as NASGRO) 
became available in 1986 from the NASA Johnson Space Center. The pro- 
gram was developed under the guidance of the NASA Fracture Control 
Analytical Methodology Panel and contains stress intensity factor solu- 
tions to a number of commonly used crack geometries. Service life 
calculations are performed with the modified Forman equation which 
reduces to the Walker or Paris equation depending on material constants 
used. 


NASA/FLAGRO is menu driven and prompts the user for information in 
a serial manner. After selecting the type of analysis desired, such as 
Safe Life, the user answers a series of questions and enters data 
depending on the particular path taken. Generally, the program operates 
serially, requiring the user to follow the same path and answer a 
number of basic questions before each execution. 

The computer program NASCRAC is being developed by Failure Analysis 
Associates under contract to NASA at Marshall Space Flight Center. For 
safe life analysis, NASCRAC has basically the same capabilities as 
NASA/FLAGRO, although implemented differently. Generally, stress inten- 
sity factors are obtained from influence function solutions to various 
geometries for which exact solutions do not exist. NASCRAC enables the 
user to select any one of several commonly used crack growth equations; 
including the Paris, Walker, and modified Forman equations. 

NASCRAC is similar to NASA/FLAGRO in that the program is menu 
driven and the user answers questions and enters data in response to 
screen prompts. With NASCRAC however, the user is not required to 
answer a series of questions before each execution. The user may 
randomly select only those menu items relating to the particular solu- 
tion desired. 
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NASA/FLAGRO SAFE LIFE FEATURES 

NASA/FLAGRO features which can affect safe life calculation: 

(1) For surface cracks with constant amplitude loading, AK is multiplied 
by a crack closure factor 

0.9 + 0.2 R^ - 0.1 R^ ; R > 0 


; R < 0 


This can increase fatigue life. 


(2) fatigue threshold is calculated using, 

AK . = (1 - C R)^ AK 
th o o 

To be conservative, let C = d = 1 for R 0 

o 


AK , = (1-R) AK 
th ^ o 


For small cracks, a ^ 0.025 in., AK^ = 0 


(3) Input of 


K - plane strain fracture toughness 

Kj^ - fracture toughness for an elliptical crack 


A^, - fit parameters 


To calculate K^ - critical stress intensity. 


a) t 


■ ■ (¥)■ 


) ■ "ic » ^ \ > 


K is incorporated into the modified Forman equation to accelerate da/dn 

c 

as K„ is approached. 
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NASCRAC SAFE LIFE FEATURES 


(1) Crack growth equations 

a) Modified Forman Analytical Comparison 

b) Walker Comparison to test data on both codes 

(2) Piecewise Linear Approximation method used. 

(3) K-solutions are based on influence functions with the default 
order of accuracy. 
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ANALYSIS NOTES 


Surface Flaws 


NASCRAC 


NASA/ FLAG RO 


Uses K, value to accelerate DA/DN per the Forman equa- 
Ic 

tion and defines failure when AK > K. , where K, is 
11 . Ic 

manually input. 

Uses K value calcualted from K, and other variables 
c Ic 

to accelerate DA/DN per the Forman equation when 

AK > K, , where K. is a material constant for surface 
.1 le’ le 


Growth Rate Equations 


NASCRAC Uses the following crack growth rate equations: Paris, 

Modified Forman, Walker, Collipriest, and Hop Rau. 

NASA/FLAGRO Primarily uses the Modified Forman Equation but the 

Paris and Walker equations could be used. 


K Values 
c 

NASCRAC K is used in the Modified Forman equation but K, is 

C -LC 

the controlling cutoff value. 

NASA/FLAGRO For ^ 0, NASA/FLAGRO uses a larger than 

for thin material (Reference 8) 

K /K, = I - B, e~” 

c Ic k 

when B, = 0, K = K, 
k c Ic 
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COMPARISON ANALYSIS CHART 


Type of Geometry 
Through Center Crack 


Through Edge Crack 


Parameters 

W = 10.0 
t - 0.25 
4130 Steel 
a^ = 0.05 

= 50 Ksi 

W = 10.0 
t = 0.25 
4130 Steel 
= 50 Ksi 

= 50 Ksi 
a^ = 0.05 


Type of Run NASGRO 

*R = 0 Tension Only X 

*R = “1 Tension Only 
Closure 
No Closure 


* R = 0 Tension Only 


*R = -1 Tension Only 
Closure 
No Closure 

*R = 0 Bending Only 

*R = -I Bending Only 
Closure 
No Closure 

*R = +0.5 Tension 
Closure 
No Closure 

*R = -0.5 Tension 
Closure 
No Closure 


NASCRAC 

X 


*See analysis results. 


XX XX XX XXX X XX 
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Type of Geometry 


Parameters 


Through Edge Crack 
(Continued) 


a = 0.25 


W = 1.75 

t = 0.A4 

D = 0.375 

B = 0.83 

4340 Steel 

0 = 59 Ksi 

t 

0 . = 37 Ksi 
b 

a = 0.25 

Through Crack at 

Pin Loaded Lug U = 5.0 

t = 0.25 
D = 0.5 
0t= 150 Ksi 
4130 Steel 
a = 0.05 
a = 0.10 
a^ = 0.25 


Through Crack at 
Pin Loaded Hole 


*See analysis results. 


Type of Run 


NASGRO 


NASCRAC 


R = +0.5 Bending 

Closure X 

No Closure X 

R = “0.5 Bending 

Closure X 

No Closure X 

R = 0 Tension X 

R = “1 Tension 

Closure X 

No Closure X 


*R = 0 Tension 
+ 

Bearing 


R = 0 
R = 0 
*R = 0 


X 

X 

X 


X 

X 

X 


XX X XX XX 
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Type of Geometry 


Parameters 


Surface Flaw Center 

Crack Specimen *Test Spec. No. 

W = 4 
t = 0.50 
= 84 Ksi 

a^ = 0.06 

a^/2c^ = 1/2 

Ti = 6AL-4V 

*Test Spec. No. 
W = 4 
t - 0.50 
= 57 Ksi 

a. = 0.06 
1 

a. /2c. = 1/2 

i' 1 

Ph-13-8M 


*See analysis results. 


Type of Run 


NASGRO 


NASCRAC 


62 


R = +0.05 


X 


X 


5576 


R = +0.05 


X 


X 



GEOMETRY 
MODELS MAacBAc 


NAagnp 





6M 

Wt* 


THROUGH 

CENTER 

CRACK 


TC01 




*0 

9l 


TTT 


M 

-w- 


i-^4. 




TC02 


THROUGH 

EDGE 

CRACK 


omlltfflS 



nnqjjniB 


202 


Ozz 

cm 






0Z2 

203 


Sg ♦P/Wt 


t t t t t 



THROUGH CRACK 
AT HOLE 
PIN LOADED 


TC03 


zz 
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ANALYSIS RESULTS 
THROUGH CENTER 
CRACK 

f w=10" 

' t=0.25" 

2 a ^0.05" 
4130 STEEL 

=50KSI 

1) R=0 

I 2) R=-1 



NASGRO MODEL TYPE TC01 NASCRAC MODEL TYPE 202 


FIRST. CA S E R=Q 


NASGRO 

K„ax=®0 ‘*0 KSiTIn @ 20,173 CYCLES 2a=1.584" 
NASCRAC 

Kmax=80 KSITTn @ 20,176 CYCLES 2a=1.60” 


SEC ON D CAS E R =- 1 

NASGRO 

CLOSURE: Kmax= 80 ‘'2 KSI flN @ 16,401 CYCLES 2a=1.584" 
NO CLOSURE: Km,A}p80.04 KSI flFl @ 4,459 CYCLES 2a=1.58" 

NASCRAC 

CLOSURE: K„ax=80 KSTTn @ 4433 CYCLES 2a=1.60" 

NO CLOSURE: SAME AS CLOSURE 
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THROUGH EDGE 
CRACK 

Si 



sic:^ 







W e 10" 
t s 0.25" 

C| B 0.05 
4130 STEEL 
O|c50KSI OR 

1) R > 0 

2) R B -1 

3) R B ^.5 

4) R B -0.5 


T 

t 


0|y B 50 KSI 

1) R B 0 

2) R B.1 


NASGRO MODEL TYPE TC02 NASCRAC MODEL TYPE 203 

TENSION ONLY 

FIRST CASE R B 0 
NASGRO 

K„ax = ® CYCLES SfB 0.610" 

NASCRAC 

Kmax ~ ® CYCLES SfB 0.655" 


SECOND CASE R b -1 



NASGRO 

CLOSURE: 

Kmax ‘ ^ cycles 

8f b0.61" 

NO CLOSURE: 

Kmax ” ^ ksi VTfT <s> 2i4s cycles 

Sf b0.61" 

NASCRAC 

CLOSURE: 

Kmax - ksi VIFT & 4439 cycles 

Si b b.655 

NO CLOSURE: 

same as closure 
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THROUGH EDGE CRACK 
TENSION CASES CONTINUED 


THIRD CASE R=.^0.5 

CLOSURE KSIVTFT @ 27,485 CYCLES a=0.61" 

NO CLOSURE Kmax= 80.05 KSIOT @ 43,772 CYCLES a=0.61'* 


NASCRAC 

CLOSURE Kj^ax= 83.22 KSlfiFT @ 43,516 CYCLES a=0.655" 
NO CLOSURE SAME AS CLOSURE 


FOURTH CASE R=-0.5 
NASGRO 

CLOSURE K„ax=® 000 KSl/IrT® 8757 CYCLES a=0.61" 

NO CLOSURE Kmax=® 0-00 KSI flFT@ 4010 CYCLES a=0.61" 


NASCRAC 

CLOSURE a=0.655” 

NO CLOSURE SAME AS CLOSURE 



THROUGH EDGE CRACK CONTINUED 
BENDING ONLY 


FIRST CASE R=0 
NASGRO 

KsiVnr @ 10,830 CYCLES a =0.73 in. 

NASCRAC 

Kmax=®‘* -363 KsiVnr @ 10,228 CYCLES a =0.721 In. 

SECOND CASE R=-1 
NASGRO 

CLOSURE: K,y,Ax=80.03 KsiVIrT @ 8846 CYCLES a =0.73 in. 

NO CLOSURE: Kmax=80.03 KsIVTiT @ 2405 CYCLES a =0.73 In. 

NASCRAC 

CLOSURE: KslVTiT @ 2271 CYCLES a =.721 In. 

NO CLOSURE: SAME AS CLOSURE 
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THROUGH CRACK 
AT PiN LOADED HOLE 


W 



W > 1.75" 

I a 0.44" 

HOLE DIAMETER a 0.375" 
EDGE DISTANCE a 0.83" 
4340 STEEL 
aTs59 KSI * 0 |,,ar 
CRACK LENGTH a 0.05" 

R a 0 


NASGRO MODEL TYPE TC03 NASCRAC MODEL TYPE 208 
NASGRO 

Kmax ■ ^ ® CYCLES 8f aO.339" 

NASCRAC 

Kmax » ksi /in eeos cycles af 30.492” 
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THROUGH CRACK 
AT LUG 



NASGRO MODEL TYPE TC04 NASCRAC MODEL TYPE 209 

WIDTH=5.0” 

THICKNESS=0.25" 

4130 STEEL 
DIAMETER OF HOLE 0.5" 
at =0.25" 

=150KSI 

NASGRO RESULTS 

Kmax=®0-71 KS\fW@ 64,426 CYCLES af=1.99” 


NASCRAC RESULTS 

Kmax=®0 ksi WFT @113,649 CYCLES a/ =2.184” 
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WALKER CONSTANTS FOR PART 
THROUGH CENTER CRACK ANALYSIS 


MATERIAL: PH 13 - 8 Mo TEST CASE NO. 5576 

c 7.63 X 10*^^ ja IN/CYCLE 

m 1.0 

n 3.54 

8 KSlV IN 
K,g 100 KSlV IN 

MATERIAL: TI - 6AL - 4V TEST CASE NO. 62 

c 2.914 X 10*^2 /A IN/CYCLE 

m 0.04435 

n 4.51 

AKth 4.5 KSlV IN 

K,g 70 KSlV IN 
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CRACK SIZE INCHES 
DEPTH a HALF LENGTH c 




a or c 


NASA/FLAGRO PART THROUGH CRACK ANALYSIS 


DATA FROM "NASGRO #5576" 



TRANSITION TO 1-D SOLUTION TC01 AT 43,531 CYCLES 

a=0.5" 1=0.5" c=0.62" 

FAILURE OCCURED AT 43,843 CYCLES K max=100.3 KSIflN c=0.798" 










DEPTH a OR HALF LENGTH e 



0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10,000 

CYCLES 



NASA/FLAGRO PART THROUGH CRACK ANALYSIS 

DATA FROM "NASGRO #62" 



2000 


4000 


6000 


CYCLES 

FAILURE OCCURED AT 7855 CYCLES asO.358** CsO.4" Kmax=85.09 KSI ^TiFT 






OBSERVATIONS 


NASCRAC and NASA/FLAGRO are both user friendly fracture mechanics 
analysis codes. Both programs offer a wide variety of crack geometries. 
Material property data can be read in from a resident file or from user 
defined input. Load spectra data for the constant amplitude loading 
cases were utilized easily in both programs. 

For the through-crack comparison analysis the Modified Forman 
equation was used and for the part-through crack analysis the Walker 
growth rate equation was used. 

For the through-crack analysis with an R ratio of zero, results 
showed good correlation between the two codes except for the through- 
crack at a lug solution. For R = -1, -K).5, -0.5, NASA/FLAGRO calculates 
an m value that is not readily known to the user; it must be hand cal- 
culated for use in NASCRAC, By specifying the non-closure option, m is 
automatically set to zero. The non-closure option gave the most con- 
servative results in NASA/FLAGRO. For R = -1, 40.5, -0.5, changing the 
m value in NASCRAC had no effect on the results, the m value has been 
permanently set to some prescribed value. The NASCRAC results for the 
through-crack analysis for R = -1, +0.5, -0.5 were in the range of the 
NASA/FLAGRO results for the non-closure option. 

For the part through center crack analysis, both programs gave com- 
parable results, particularly with specimen No. 15576 where the crack 
grew through before failing, but both programs showed failure before 
breakthrough for specimen No. 62 which was different from the results 
of the test. 

The comparison analysis between the two programs is an on-going 
effort for our analysis team. Other types of solution methods and 
problems are scheduled to be studied in the future. 
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